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Abstract 

This paper studies an infinite-server queue in a random environment, 
meaning that the arrival rate, the service requirements and the server 
work rate are modulated by a general cadlag stochastic background pro¬ 
cess. To prove a large deviations principle, the concept of attainable pa¬ 
rameters is introduced. Scaling both the arrival rates and the background 
process, a large deviations principle for the number of jobs in the system 
is derived using attainable parameters. Finally, some known results about 
Markov-modulated infinite-server queues are generalized and new results 
for several background processes and scalings are established in examples. 
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1 Introduction 

The infinite-server queue is one of the fundamental models in queueing theory. 
Its distinguishing feature is the presence of an infinite number of servers, so that 
jobs are served independently and there are no waiting times. This leads to 
explicit formulas for many quantities of interest, especially for M/M/oo queues, 
where jobs arrive according to a Poisson process and the service requirements 
have an exponential distribution. In practice, however, one often observes time- 
varying arrival intensities, service requirement distributions and server work 
rates. This calls for adequate modeling. 

A natural way to incorporate time-dependence is to consider an M/M/oo 
queue in a random environment. In this case there is an independent background 
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process that modulates the arrival rate, the service requirement distribution and 
the work rate of the servers. 

Model. In this paper, we study the case where the background process is a 
general stochastic process J whose paths are right-continuous and have finite 
left limits, i.e. J has cadlag paths. The process J modulates the arrival rate, the 
service requirement distribution and the server work rate in the following way. 
When J is in state cc, jobs arrive according to a Poisson process with intensity 
A(*). Upon arrival, a job draws an independent service requirement from an 
exponential distribution with parameter k{x) if J is in state x when the job 
arrives. Then the service requirement of the job is processed by a server, whose 
work rate is n(x) while J is in state x. Immediately after its service requirement 
has been processed, a job leaves the system. 

Main result. The main result of this paper is a full large deviations principle 
(LDP) for the transient number of jobs in the system, under a scaling of the 
arrival rate and the background process. To arrive at this result, we first show 
that the number of jobs in the system at time t > 0 has a Poisson distribution 
with random parameter </> t (J). Then we scale A H > nX and we scale J i—► J n 
such that the normalized random parameter 4>t{J n ) satisfies an LDP. Under this 
scaling, we derive the LDP for the transient number of jobs in the system. 

Literature. The amount of literature on infinite-server queues in a random 
environment is quite small. Moreover, almost all papers on this topic (with 
notable exception 0) study Model I or Model II (cf. [5]). In both models, jobs 
arrive according to a Poisson process with intensity A(x) when the background 
process is in state x. In Model I, service requirements have a standard exponen¬ 
tial distribution and servers work at rate n(x) when the background process is 
in state x. This is equivalent to the jobs being subject to a modulated hazard 
rate. In Model II, service requirements have an exponential distribution with 
parameter n(x) when the background process is in state x and servers work at 
constant rate 1. 

An early reference is m, which analyzes Model I when the background 
process is a continuous-time Markov chain. Important results in El are a 
recursion for the factorial moments of the number of jobs and the observation 
that the steady-state distribution is not of some ‘matrix-Poisson’ type. 

Other important results can be found in [7], which studies Model I when 
the background process is a semi-Markov process with finite state space. The 
crucial observation in [7 is that the stationary number of jobs has a Poisson 
distribution with a random parameter that is determined by the background 
process. Moreover, the factorial moments of the number of jobs are computed 
via a recursion. These results are generalized in |12j . 

The observation in [7] is used to obtain time-scaling results in both the 
central limit regime and the large deviations regime. In the central limit regime, 
[2] and [3] derive central limit theorems for Markov-modulated infinite-server 
queues for several models and scalings. In this regime, the so-called deviation 
matrix (cf. i6]) plays an important role. In the large deviations regime, [3] and 
[5] compute optimal paths to obtain rate functions under a linear scaling of the 
arrival rates, given that the background process is an irreducible continuous-time 
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Markov chain. The former studies Model I, whereas the latter studies Model II 
for a class of service requirement distributions that includes the exponential 
distribution. 

As mentioned, we show that the number of jobs in the system has a Poisson 
distribution with a random parameter, which can be interpreted as a mixture 
of Poisson distributions. In [T], an LDP is derived for mixtures that satisfy 
certain assumptions. However, apart from the assumption that the normalized 
random parameter satisfies an LDP, these assumptions are either superfluous or 
too restrictive in our case. In particular, we do not assume that the sequence of 
measures induced by the normalized random parameter is exponentially tight, 
so we cannot use the arguments in pQ. Hence, we need a different approach to 
obtain an LDP. 

Contributions. In more detail, the contributions of this paper are the follow¬ 
ing. We generalize known models by considering a general cadlag background 
process instead of a semi-Markov background process with finite state space. 
Moreover, in our model the background process modulates both the service 
requirement distributions and the server work rate, whereas previous papers 
considered models in which either the service requirement distributions or the 
server work rate was modulated. In particular, our model generalizes Model I 
and Model II. 

Using elementary arguments, we show that in this model the transient num¬ 
ber of jobs has a Poisson distribution with random parameter. We scale the 
arrival rate linearly and we scale the background process such that the nor¬ 
malized random parameter satisfies an LDP. Under this scaling, we obtain a 
full LDP for the number of jobs in the system. To the best of our knowledge, 
this is the first time that a full LDP is presented for modulated infinite-server 
queues. To prove the LDP, we introduce the concept of attainable parameters 
and use a variation on Varadhan’s Lemma. These tools enable us to avoid the 
assumptions in pQ. 

The theory is illustrated by examples that show rate functions that cannot 
be obtained via background processes with finite state space. Additionally, we 
show that completely different background processes may lead to the same LDP, 
even in highly nontrivial cases. 

Organization. The rest of this paper is organized as follows. In Section 2, we 
describe the model and provide some of its basic properties. Additionally, we fix 
some notation. In Section 3, we introduce the concept of attainable parameters 
and prove an LDP for the number of jobs in the system. In Section 4, we 
show that the rate function corresponding to this LDP has a simple description 
when we do not scale the background process. As an illustration, we work 
out some examples. In Section 5, we give examples in which we do scale the 
background process. In Section 6, we briefly discuss the results and point out 
some topics for future research. The appendices provide some technical details 
about the number of jobs in the system (Section A), continuity in Skorokhod 
space (Section B) and properties of Poisson random variables (Section C). 
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2 Model and problem description 

We study an infinite-server queue with modulated arrival rates, service require¬ 
ments and server work rates. The precise mathematical setup of the model and 
some of its basic properties are provided in Section [A] In words, the model may 
be described as follows. 

Let (J(t)) t>0 be a cadlag stochastic process with state space £, which is 
assumed to be a metric space. We will refer to the process J as the background 
process or modulating process. While the background process is in state x E £, 
jobs enter the system following a Poisson process with intensity A(x) > 0. 

When job k enters the system, it draws a service requirement from an inde¬ 
pendent exponential distribution with parameter re(y) if the background process 
is in state y € £ upon its arrival. Server k processes this service requirement 
at rate p(z) while the background process is in state z € £. Job k leaves the 
system when its service requirement has been processed. 

We denote a modulated infinite-server queue by the quadruple (J, A, k,/ z). 
Additionally, we denote the number of jobs in this system at time t by M(t). 
In Section El it is shown that M(t) has a Poisson distribution with random 
parameter 


pt 

MJ)= / A(J(s))e- K(J(s)) ^> (JW)dr ds. (1) 

Jo 

This will turn out to be a crucial property in this paper. 

We are interested in events with an unusual number of jobs in the system. 
More precisely, we would like to prove an LDP for the number of jobs in the 
system. A sequence of probability measures {r n } ngN is said to satisfy an LDP 
with rate function p if there exists a lower semi-continuous function p: X —> 
[0, oo] such that 


limsup — logr„(F) < — inf p(a) 

n—>oo Tl a£.F 


for all closed sets F and 


liminf — log r n (G) > — inf p(a ) 

n—> oo 77, a£G 

for all open sets G, where each t„ is defined on the Borel er-algebra of the 
topological space X. A sequence of random variables is said to satisfy an LDP 
with rate function p if the sequence of measures induced by the random variables 
satisfies an LDP with rate function p. Importantly, we do not assume that p is 
a good rate function, i.e., we do not assume that p has compact level sets. 

As mentioned, we would like to prove an LDP for the number of jobs in 
the system. To analyze this problem, we will scale the arrival rates via A(ir) H>• 
nX(x), i.e., we linearly speed up the arrivals. In addition, we will scale the 
background process via J i-> J„. Formally, scaling A(a;) H > nX(x) and J i —> J n 
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means that we start with an infinite-server queue (J, A, n, p) and then consider 
the sequence of infinite-server queues {( J n , nX , k. /r)} ngN - 

Given the scalings X(x) K > nX(x) and J i—>■ J ni we denote the corresponding 
number of jobs in the system by M n (t). It follows immediately from equation (fTj) 
that M n (t) has a Poisson distribution with random parameter 



The normalized random parameter (f>t{J n ) induces a sequence of probability 
measures {i/„} ngN on v ^ a v n(B) = V((j> t {J n ) £ B) for Borel sets Bel. 

We will assume that the sequence of probability measures {^n} ngN satisfies 
an LDP with rate function ip. Note that {z 7 i}„ 6 n trivially satisfies an LDP 
when v n = tVi+i for all n £ N, so this assumption covers the case in which the 
background process is not scaled. 

Given the scaling, we denote the number of jobs in the system at time t by 
M n (t) and consider the normalized random variable i M n (t ). Our goal is to 
prove an LDP for -M n (t) and to describe the corresponding rate function. 

Throughout this paper, we will also use the following notation. We denote 
the closure of a set A by cL4. We write B(x,e) for the open ball with center 
and radius e > 0 and B[x, e] for its closure. For notational convenience, 
we will sometimes write R + for [0, oo), B + (x,e) for B(x,e) D R+ and B + [x,e] 
for B[x, e] D R+. As is customary, we define exp(—oo) = 0 and log(0) = —oo. 

3 A large deviations principle 

In this section we will prove an LDP for the number of jobs in the system under 
a scaling of the arrival rates and the background process, i.e., we will prove an 
LDP for i M n (t ). It will turn out that so-called attainable parameters determine 
the rate function corresponding to the LDP. 

Definition 3.1. Given a scaling J i —> J n , a real number 7 £ [0,oo) is called an 
attainable parameter at time t > 0 if for all e > 0 there exists N e £ N such that 
G B{ 7 ,e)) = v n {B> 0 for all n > N t . The set of all attainable 
parameters at time t is denoted by lZ(t). 

The intuition behind attainable parameters is as follows. The number of 
jobs in the system has a Poisson distribution with a random parameter that is 
completely determined by the background process. Basically, the background 
process samples the Poisson parameter. A real number 7 is an attainable param¬ 
eter if, for all n large enough, the scaled background process samples parameters 
close to 7 with positive probability. 

As mentioned before, we will prove an LDP for —M n (t) by scaling X(x) *-> 
nX(x) and J 1 —>■ J n such that the sequence of probability measures {iz ra } ngN sat¬ 
isfies an LDP with rate function ip. The rate function I: R —► [0, 00 ] governing 
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the LDP for i M n (t) is given by 

I ( a ) = inf [£(r,a) + ip(j)}, ( 2 ) 

76 n(t) 

where £( 7 ; •) is the Fenchel-Legendre transform of the Poisson cumulant gener¬ 
ating function with parameter 7 . It will turn out (cf. Lemma l3.2|) that 

1(a) = inf [£( 7 ; a) + V >( 7 )] = inf _ [^(75 a) + ip(l)\. (3) 

767 <,(() 76 {y<oo} 

However, we will take the infimum over lZ(t) rather than over {ip < 00} to stress 
that attainability of parameters is the crucial property for proving the LDP. 

Before we can give the proof, we have to settle some technical details. First, 
it is not immediately clear whether the function I is indeed a rate function or 
even whether I is well defined. In particular, it is not clear whether lZ(t) is a 
non-empty set. However, the assumption that the sequence { 77 i} ngN satisfies 
an LDP implies that 7 Z(t) is non-empty, as the following lemma shows. 

Lemma 3.2. Let the scaling J 1 —>• J n he such that {Pn} ngN satisfies an LDP 
with rate function ip. Then 7 Z(t) is a non-empty closed subset of [ 0 , 00) and 
{ip < 00} C 7 Z(t). 

Proof. Suppose that 7 £ R \ 7 Z(t). Then there exists e > 0 such that for all 
n £ N there exists k n £ N such that k n >n and Vk n (B( 7 ,e)) = 0. This implies 
that 77 ( 7 , e) C R \ 7 Z(t), so 7 Z(t) is closed. Moreover, we must have 

lim inf — log e)) = — 00 = — inf ip(a), 

n—too n aeB(j,e) 

so ip (a) = 00 for all a £ B{j, e). Then R \ 7 Z(t) C {ip = 00} and {ip < 00} C 
lZ(t). The fact that ip is a rate function implies that {ip < 00} is non-empty. 
The statement of the lemma follows immediately. □ 

From the previous lemma it follows that I is a well defined function. The 
fact that J is a rate function is implied by Proposition 1C.51 and the functions l 
and ip being rate functions. 

The next lemma is a variation on Varadhan’s Lemma. Contrary to Varad- 
han’s Lemma, it does not require that a given function / is continuous. Instead, 
it requires that a weaker condition is fulfilled. We will use this lemma to obtain 
the large deviations upper bound, by applying it to functions / of the form 
described in Proposition 1C .41 

Lemma 3.3. Let X be a topological space and let {£n} ng ^ be a sequence of 
measures defined on its Borel o-algebra. Suppose that {£n}„ g N satisfies an LDP 
with rate function g. Let /: df —>■ [— 00, 0 ] be a Borel measurable function such 
that /^ 1 ([a, 6]) is a closed set for all a, 6 £ (— 00, 0 ] satisfying 

sup [/(a;) — g(x)} < a < b < 0 . 

xGX 


6 




Then it holds that 


limsup —log f e n f( x ' > £ n (dx) < sup[/(a;) — p(x)]. 
ra—> oo n Jx xeX 

Proof. This follows immediately from HU Lem. 2.2] □ 

With these technical details settled, we can prove the following LDP. Its 
proof exploits two elementary observations. First, conditional on a value of the 
random parameter the number of jobs in the system has the same distri¬ 

bution as the number of jobs in the system in the M/M/oo setting. Second, the 
number of jobs in the system in the M/M/oo setting has the same distribution 
as a sum of i.i.d. Poisson random variables. Combined with some analytical 
results, these observations enable us to prove the LDP. 

Theorem 3.4. Consider a modulated infinite-server queue (</, A, k,/i) as de¬ 
scribed in Section 0 Scale X(x) K > n\(x) and J K > J n such that {^n} ng is} 
satisfies an LDP with rate function ip. Then the rescaled number of jobs in the 
system —M n (t) satisfies an LDP with rate function L as defined in equation ([2)1 . 
so 

limsup — logP[ —M n (t) G F ) < — inf 1(a) (4) 

n->oo U \H ) a£F 

for any closed set Fcl and 

liminf — logPf —M n (t) G G ) > — inf 1(a) (5) 

n-¥ oo n \n ) aeG 

for any open set Gel. 

Proof. For A > 0, let Po(X), Pi (A), P 2 (X),... denote a sequence of i.i.d. random 
variables that have a Poisson distribution with parameter A. Let F C 1 be a 
closed set and let G C R be an open set. 

To prove the upper bound (HI), recall that M n (t) has a Poisson distribution 
with random parameter ncpt(Jn)- Then we may write 

limsup — logPf —M n (t) € F ) = limsup —log / P( — Po(n^) £ F] v n ( dy) 
n—too Tl \Tl J n—too Tl J[0,oo ) / 

= limsup -log [ P( - V P* ( 7 ) G F ) v n (dry) 
n^oo n 7[0,oo) \ n J 

< limsup - log [ 2e n[ - ini ^ Fe{T ’ a)] i/ n (d 7 ) 
n-¥ 00 n J [0,oo) 

= limsup - log f e n[-i»f. g F^(7ia)] ^(dy). 
n->oo Tl J [0,00) 

The inequality above is an immediate result of the proof of Cramer’s Theorem 
in K as provided in 10). 
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According to Proposition 1C.41 the function 7 i-»- — inf ae p ^( 7 ; a) satisfies 
the assumptions of Lemma 13.31 Moreover, {^ n } ngN satisfies an LDP both in R 
and in [0,oo) with rate function 'ip (cf. [TO,, Lem. 4.1.5]). Hence, we may apply 
Lemma 13.31 to obtain 


lim sup — log P 

n—too Tl 


-M n (t ) e F 


< lim sup — log 

n—^oo Tl 


a n[- inf aeF e(T,a)} 


[0,oo) 


e ' 1 i'n(d'y) 


< sup 

7E[0,oo) 


- inf £( 7 ; a) - ip( 7 ) 

aGF 


= -inf inf [t{r, a) + ^( 7 )] 

a(zF 7G [0,oo) 

= - inf inf [£( 7 ; a) + ^( 7 )] 

Q6F 7G -R(t) W/J 

= — inf 1(a). 

a£F 


The fact that we only have to consider the infimum over lZ(t) follows from 
Lemma 13.21 This proves the upper bound. 

To prove the lower bound ©, let A £ 7 Z(t) and e > 0. Define A e = 
max{0, A — e} and A+ = A + e. By definition of the set 7 Z(t) there exists N e such 
that P (4>t(Jn) £ B( A, e)) > 0 for all n > N e . 

Fix x £ G. Because G is open, there exists S > 0 such that B(x,5) C G. 
Observe that 


—M n (t) £ G ] > 

n 1 




—M n (t) £ B(x,S) 
n 


M n (t) £ B(x,6 ); (p t (J n ) £ B( A, e) 


M n (t) £ B(x,6) 


4>t(Jn) £ B(X,e) )P (4>t(J n ) £ B( A,e)) 


for all n > N e , where the equality follows from the fact that P ((f> t (J n ) £ B(X,e)) > 
0 for all n> N e . Then we get 


1 


n—too n 


1 


lim inf — logP —M n (t) £ G ] > 


n—> 00 fi 


lim inf — log P ( — M n (t ) £ B(x,5) 


(pt(Jn) £ B(X, e) 


+ liminf-logP(0 t (J n ) £ B( A, e)). 

n—too Tl 

Recall that (pt(J n ) satisfies an LDP with rate function ip, so 


lirninf — logP (cpt(Jn) £ B( A, e)) > — inf ip (a) 

n-¥ 00 n a£B(\,e) 
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by assumption. Moreover, it holds that 


1 


n—>oo 77 , 


liminf — logP —M n (t) G B(x, S) 


n 


n—> oo 77 , 


liminf — logP —M n (t) £ B{x , 8) 


lim inf inf — 

n-> oo 7 eB(A,e)nR + 71 


<Pt(Jn) € B(X,e)j — 
4>t(Jn) € B( A, e) n > 


i= 1 


- inf ^(7; a) 

aGB(x,(5) 


7 e{Ar,A+} 

The equality above is established in Proposition [(13] Combining the results, we 


obtain that 


— M n (t) £ G ) > min 

' 76{A7,A+} 


- in f . At, a) 

a£B(x,S) 


— inf ip (a). 

o 6 fl(A,e) 


This holds for all e > 0 and small enough 8 > 0. Taking limits, we get 


lim 


— inf f(y; a) 


40 7 e{Ar,At}L a£B(x,S) 
thanks to Proposition [C3] and 


= — inf £(A; a) 

a(zB(x,5 ) 


fim inf ip(a) = ip( A), 

because ip is lower semi-continuous. Similarly, we get lim^o inf a eB(a;,d') ^(A; a) = 
^(A;x). Hence, it follows that 


1 


) > lim lim 

min 

/ <54.0 ej.0 

_ 7 e{A e -,A+} 


— inf £( 7 ; a) 

a(zB(x,5) 


— inf ip (a) 

a£B( A,e) 


= — [^(A; a;) + V>(A)]. 

Since x £ G and A £ TZ(t) were arbitrary, we obtain 

—M n (t) £ G ) > sup sup [-[^(A;cc) + -0(A)]] 


n 


cl^lG A£7?.(i) 

= — inf 1(a ), 
aeG v ’ 


which completes the proof. 


□ 


The proof of Theorem l.T 41 contains familiar elements. First, the upper bound 
is proved using a Chernoff bound combined with a variation on Varadhan’s 
Lemma. Second, the lower bound is proved by considering ‘the most likely 
of all unlikely scenarios’, which is similar to the method used in [3j and [5]. 
However, the proofs there relied on properties of irreducible continuous-time 
Markov chains and the computation of optimal paths, whereas we consider 
general cadlag background processes via attainable parameters. 
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4 Examples: unsealed background processes 


Given the scaling A H > nX and J i-> J„, Theorem 13.41 provides a full LDP 
for — M n (t ) and describes the corresponding rate function. In the upcoming 
examples we will consider cases in which the background process is not scaled 
and we will use Theorem m to verify or extend known results and to obtain 
new results. 

Throughout this section we will assume that the background process is not 
scaled, i.e., J n = J for all n G N for some cadlag stochastic process J. The 
following lemma is trivial, but plays a central role in this section. 

Lemma 4.1. If J n — J for all n G N, then the sequence {( Pt(Jn)} n satisfies an 
LDP with rate function ip. In this case it holds thatlZft) = {if < oo} = {ip = 0}. 

Hence, when the background process is not scaled, we have the special prop¬ 
erty that 7 Z(t) = {ip = 0}. This will enable us to compute explicit rate functions 
in the examples. In these computations, we will extensively use the following 
properties of the rate function I and properties of step functions in Skorokhod 
space. 

Recall that the rate function I is given by 

1(a) = inf ' Mr, a) + iph)], 

7 eiz(t) 


and that lZ(t) = {ip = 0} (see Lemma r4.il) . Hence, we get 

/(«)= inf %;a). 

7 en(t) 


( 6 ) 


In this case, we can give a simpler and more explicit description of I, using the 
following properties of the function I. 

For 7 > 0, the function ^( 7 ;-) is the Fenchel-Legendre transform of the 
Poisson cumulant generating function with parameter 7 and is given by 


00 a < 0; 

?(r,a)={ 7 a = 0; 

7 — a + a log(a/ 7 ) a > 0 . 


(7) 


For 7 = 0 and a > 0 , we understand that 7 —a + alog(a/ 7 ) = 00. An important 
observation is that the following inequalities hold for 0 < 71 < 72 < 00: 


a) < ^( 72 ; a) 
^(7H o) > e(j 2 ;a) 


Va G [0,71]; 
Va G [72, 00). 


( 8 ) 

(9) 


See Figure Q] for an illustration. 

Because in the present case I is just an infimum of Poisson rate functions, 
these inequalities imply that / has some special properties. They are described 
in the following proposition. 
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7i 72 


Figure 1: Graphs of the functions •) and ^( 72 ; ■) for 0 < 71 < 72 < 00 



Figure 2: Visualization of the function / in Example 14.31 


Proposition 4.2. In the present case, I {a) =0 if and only if a £ IZ(t). If 
1(a) > 0 for some a £ K., then exactly one of the following three scenerios is 
true: 

1. a < C- = milZft) and 1(b) = £(c_; b) for all b £ (— 00 , c_]; 

2. a > c+ = swplZ(t) and 1(b) = £(c+;&) for all b £ [c+, 00 ); 

3. the previous two cases do not hold and 1(b) = min{T(c_; 6 ),f(c+; b)} for all 
b £ [c_,c+], where C- = sup(7 Z(t) D (—00 ,a)) andc+ = sup(7?.(t) D (a, 00 )). 

Proof. It follows immediately from equations ([Gj and 0 that 1(a) = 0 if and 
only if a £ lZ(t). Hence, 1(a) > 0 implies that the distance of a to lZ(t) is 
strictly positive, since 7 Z(t) is closed. The three scenarios now follow from the 
inequalities (0 and 0 . □ 

The previous proposition may seem rather abstract. To get some intuition, 
the following example describes a typical rate function. 

Example 4.3. Suppose that 7 Z(t) = [a,j3\ U [ 7 , <5] for some 0 <a</?< 7 < 

S < 00 . Then the function I looks like the graph shown in Figure [2] it equals 0 
on the intervals [a,/3] and [7,(5], whereas it equals the minimum of £(/3; •) and 
£( 7 ; •) on the interval (/3, 7 ) in between. On the interval (— 00 , a] the function I 
equals £(a; ■) and on the interval [(5, 00 ) the function / equals £(S; •). 
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To compute 7Z(t), it is often convenient to use the following properties of step 
functions in D([0, oo); £). (For the definition of a step function, see Section iBl) 
Recall that the set of all step functions in D([0, oo); £) is denoted by <S([0, oo); £). 

Lemma 4.4. If £ <S([0, oo); £)} C Pit), then 


1 ^(i) = d{Mf)\f € 5([0,oo);f)} = e Z>([0,oo);£)}. 

Proof. This follows from Lemma 13.21 Corollary IB. 31 and Lemma IB. 41 □ 

Lemma 4.5. If {</>t{f)\f £ <S([0, oo); £)} c Hit), thenlZ(t) is a closed interval. 

Proof. It suffices to show that 7 Z(t) is convex. Let f}. , ff £ S([0,oo);£). We 
may assume that dt(/c) < <Pt (fc) ■ For x £ [0,7] we define the function g x via 

9x(s) = l{s<x}/c ( s ) T l{s>a;}/c ( s ) 

for s £ [0,oo). Clearly, g x £ 5([0, oo);£) and 

M9x) = Mfc) + (Mfc) - ^(/ca¬ 


using the continuity of the integral and applying the Intermediate Value The¬ 
orem, it follows that 


[Mfc), Mfc)] c {M9x)\x £ [0,7]} C 7^(7). 


Combined with Lemma 14.41 this implies the statement of the lemma. 


□ 


Let f c £ <S([0,oo);£) be a step function. Clearly, f c has a unique minimal 
representation {(7^, aj)}f =0 , where k £ N, 0 = to < 7i < ... < tk < oo and 
oo,..., at £ £ are such that / c (7) = for 7 £ [U, 7i+i) and i = 0,..., k — 1 
and / c (7) = for 7 £ [7fc,oo). Given this minimal representation, we define its 
truncated minimal step size by 


A fc = 1 A min {7* - 7*_i}. 
i=l,...,k 

Additionally, we define t^+i =tkVt. The truncated minimal step size and tk+i 
will be used for computing attainable parameters. 

In the upcoming examples, we would like to compute rate functions via at¬ 
tainable parameters. To compute attainable parameters, we use the following 
strategy. We fix a certain path /, often a step function. This gives us a param¬ 
eter value </>*(/)■ Then we would like to show that, with positive probability, 
the background process stays ‘close’ to /, which will imply that 4>t{f) is an 
attainable parameter. 

Staying ‘close’ to / depends on properties of £ and the background process. 
In most cases, the background process needs a little bit of room (both in time 
and in space) to jump near a discontinuity of /. This is where the truncated 
minimal step size comes in: it is an upper bound on the time we give the 
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background process for jumping near a discontinuity of a step function. The 
precise meaning of this will become clearer in the examples. 

The first example treats the familiar case of a Markov-modulated infinite- 
server queue, i.e., the case in which the background process is an irreducible 
Markov chain. This case is partly studied in [3] (Model I) and [5] (Model II). 
In the example, we recover [3j Th. 2] and [5] Th. 1]. Additionally, we generalize 
these results to our model and extend them to a full LDP. 

Example 4.6. Let J be an irreducible, continuous-time Markov process with 
finite state space £ = {l,...,d} and consider the modulated infinite-server 
queue (J, A, k,//). Given the scaling A H > nX, Theorem 13.41 (combined with 
Lemma ED shows that —M n (t) satisfies an LDP with rate function I. This 
rate function may be computed as follows. 

Note that Z?([0, oo); £ ) = <S([0, oo); £), since £ is finite. Fix any function g £ 
.D([0, oo);£) with minimal representation {{ti, aii)}i=o and take any e £ (0,1). 
Define W(g;e) as the set of all / £ Z7([0, oo); £) such that 

/ (t) = on-\ Mt £ [fi_i + f - fjrAg) Vi S {1,..., k}, 

f(t) = a k Mt £ [t k ,t k +i\. 

Now note that 

sup </> t (/) < 4>t(g) + e max A (j) 

/EWfee) 

and 


, ^ Mf) > Mg) - e . max A (j), 

f£W(g;e) 

so we can get both the supremum and the infimum arbitrarily close to by 

taking e small enough. 

Observe that P(J £ W(g;e)) > 0, thanks to the irreducibility of J. Conse¬ 
quently, 7 Z(t) = {4>t(g) | g £ B([0,oo);£)j. Then Lemma ES] implies that 7 Z(t) 
is a closed interval. Using that £ is finite, we immediately get 


Hit) = [o_,a+], 


where 0 < a_ < a+ < oo with a_ = inf geD([0jOo);f} 4> t {g) and a+ = sup geD([0 oo);£) 4> t {g). 
Now applying Proposition 14.21 it follows that the rate function / is given by 



oo 

a £ (—oo, 0); 

£(a_; a) 

a £ [0, a_]; 

0 

a £ [a_, a+]; 

£{a+;a) 

a £ [a+, oo). 


( 10 ) 


The result of the previous example depends neither on the initial distribution 
nor on the transition rate matrix of the irreducible Markov chain. Moreover, 


13 


the analysis in the previous example implies the following lemma. It shows that 
we always obtain a good rate function when the background process has a finite 
state space. 

Lemma 4.7. Let jW be a background process with finite state space £ and let 
J 1 ' 2 ' 1 be an irreducible Markov chain with the same state space. Consider the two 
modulated infinite-server queues /i) and (j^ 2 ), A, k, p). Scaling A i— > 

n\, we obtain in both cases an LDP for the number of jobs in the system with 
corresponding rate functions I ^ and I^ 2 \ Then it holds that I^\a) > I^ia) 
for all a £ R. In particular, both I^ 1 ) and fi 2 ' 1 are good rate functions. 

In the next example we will modulate an infinite-server queue by another 
Markov-modulated infinite-server queue. This setup differs from the setup con¬ 
sidered in [3] and j5J. In particular, the state space of the background process 
is countably infinite, so that we may obtain a rate function that is not good. 

Example 4.8. Consider a Markov-modulated infinite-server queue as described 
in El, i.e, a Markov-modulated infinite-server queue under the assumptions 
of Model I. Assume that neither the arrival rates nor the server work rates are 
identically equal to 0 and that the system starts empty. Let J( t ) be the number 
of jobs in this Markov-modulated infinite-server queue at time t > 0. Then J is 
a cadlag stochastic process and its state space is £ = Z^ 0 - 

Consider the modulated infinite-server queue (J, A, k. p) and impose the scal¬ 
ing A i—^ nX. Then A M n (t ) satisfies an LDP with rate function /, according to 
Theorem 13.41 and Lemma T4. II This rate function may be computed as follows. 

Recall that J stays in state m £ £ during [t, t + At] with positive probability 
for arbitrarily large At. Moreover, because neither the arrival rates nor the 
server work rates are identically equal to 0, the process J also has the following 
property. If J(t) = mi at time t > 0, then it jumps to state m 2 £ £ during 
[t, t + At] with positive probability for arbitrarily small At. 

Roughly speaking, these two properties mean that the background process 
is irreducible, in the sense that it can jump to or stay in any state during 
any time interval we would like. Of course, this is very similar to the Markov 
chain being irreducible in the previous example. Consequently, our strategy for 
determining the attainable parameters will be very similar, although there are 
some subtleties related to the state space being infinite. 

Fix any g £ <S([0,oo);£) with minimal representation {(i,,aj)} i=0 and take 
any e £ (0,1). Let W(g; e) denote the set of all / £ D([0, 00 ); £) with 

f(t) = 1 Vi G [ti_ 1 + fjrA g ,ti - §£A g) Vi G {1, ... ,&}, 

/(^) = &k Vi G 

and 


0 < f(t) < a 0 

Vt £ [0, 2^A g ), 

Oii- i A oti< fit) < ati-i V Oii 

Vf £ [ti - f£A g,ti I 


Mi £ {0,..., k — 1}, 

ak- 1 A ak < fit) < ak- 1 V ak 

Mt £ [tfc — f jfAg, tfc] 
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Then we have 


and 


sup < 4>t{g) + e max max A (j) 

feW(g;e) ie{0,...,k} je{0,...,an} 


. m f <M/) > Mff) 

feW(g;e) 


— e max max A (j). 


The two properties of the background process described above imply that 
P(J £ W(g;e)) > 0. It follows that {</>*(<?) | g S <S([0, oo); £)} C lZ(t). Write 
a- = inf s£ £)([o l00 ) ; £) <l>t(g) and a+ = sup seD([0oo) . £) </> t (g). Lemma [474] and 
Lemma m imply that lZ(t) = [a_,a+] if a+ < oo and lZ(t) = [a_,oo) if 
a+ = oo. Hence, 


if a+ < oo and 


1(a) 


oo 

a £ (—oo, 0) 

(■(a-\ a) 

a £ [0, a_]; 

0 

a £ [ft—, g+] ■ 

£(a+;a) 

a £ [a_|_, oo) 


1(a) 


{ oo a £ (—oo, 0); 

£(a_; a) a £ [0, a_]; 

0 a € [a_, oo) 


( 11 ) 


( 12 ) 


if a+ = oo. Note that / is not a good rate funtion if a+ = oo. 

The previous example only depends on the state space being countable and 
discrete and on the background process being irreducible in the sense described 
above. Consequently, the same result holds for irreducible Markov processes 
with a countable, discrete state space. 

In the last example of this section we compare rate functions that are ob¬ 
tained using two different background processes. One background process is a 
Markov chain, whereas the other background process is a reflected Brownian 
motion, which has an uncountable state space. It turns out that both back¬ 
ground processes lead to the same LDP, even though the background processes 
are completely different. Apparently, two very different modulating processes 
may lead to the same rate function for the LDP, even if the arrival rates, service 
requirements and server work rates are nontrivial. 

Example 4.9. Let £ = [0,1] be equipped with the Euclidean metric. Define 
A: [0,1] —> [0,1] by A(a:) = x, k: [0,1] —¥ [0,1] by k(x) = 1 and g,: [0,1] — > [0,1] 
by fi(x) = 1 — x. 

Let J MC be an irreducible, continuous-time Markov chain with state space 
{0,1}. Let J rBM be a reflected Brownian motion with reflecting barriers 0 and 
1. For simplicity, assume that J rBM starts in xo £ (0,1), so 

J lBM (t) =x 0 + W(t) + L(t) - U(t ) 
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for some standard Brownian motion W, lower-regulator process L and upper- 
regulator process U (see for instance 0). 

Consider the two modulated infinite-server queues (J MC , A, k, /i) and (J rBM , A, k, ji) . 
Under the scaling A H > n\, both and 1 c (t) satisfy an LDP with 

the same good rate function /, which is given by 

{ oo a £ (— 00 , 0); 

0 a £ [0, t]; (13) 

i(t; a) a £ [t, 00 ). 

The rate function for the LDP corresponding to is derived in 

Example 14.61 It is easy to see that the rate function has the form claimed 
above. 

We will show that ^M)j BM (t) satisfies an LDP with the same rate function. 

Fix g £ 5([0, 00 );£) with minimal representation {(U,ai)}*L 0 and take any 
e > 0. Define W(g; e) as the set of all / £ D([0, 00 ); £) such that 

\/t £ [ii_i + fjAg.ti - f^A s ) Vi e {1,..., k}. 

Then we get 


sup 

few(g-,e ) 


Mf) < <M3) + et + e 


and 


. inf ,Mf) > <!>t{g) - et - e. 

f£W(g-,e) 

Now observe that 

P( J rBM £ W(g-, e)) > P(a; 0 + W £ W{g; e)) > 0, 

due to the definition of J rBM and W being a Brownian motion. 

It follows that {(fit(g) | g £ <S([0,oo);£)} C 7?. rBM (t), so lZ rBM (t) = [0, i] and 
the corresponding rate function is given by the function I above. 

In this section we considered examples in which the background process 
was not scaled. As shown, this implies some special properties, which we can 
use to explicitly compute rate functions. In the next section, we will scale the 
background process, too. Although explicit computations are not possible in 
general, there are still cases for which we may derive rate functions. 


5 Examples: scaled background processes 

In this section we will give two examples in which the background process is 
scaled. In the first example, we will consider the Markov-modulated infinite- 
server queue and derive an explicit rate function under a superlinear time- 
scaling. In the second example, we will consider a new model in which the 
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background process is a Brownian motion. In this case, the rate function will 
be given as the solution of a variational problem. 

Example 5.1. Let J be an irreducible continuous-time Markov chain with 
finite state space rf} and generator matrix Q. Denote the corresponding 

stationary distribution by 7r = (7Ti,..., 7 r^). 

Consider the modulated infinite-server queue (J, A, k,/x). Define = 

Ej=1Tf/L and 



Scale A i-> nA and J i-»- J n , where J n {t) = J(n 1+e t). It is easy to see that 
scaling J i-» J n is equivalent to scaling Q i-» n 1+e Q. 

The sequence of random parameters satisfies an LDP with rate 

function ip, where 


V’O) 


0 a = Q t ] 
00 s / Q t . 


Indeed, this follows from the fact that 

liminf - logP(</> t (J„) G B(p t ,t])) = 0 

n—>-oo Ti 


and 

limsup - logP(0t(J n ) ^ B(p t ,r])) = -00 

n—>-oo ^ 

for all 77 > 0. These equalities are an immediate result from the proof of [31 
Th. 3], 

Given this LDP for {<pt{Jn)} nG ^, Theorem 13.41 implies that satisfies 

an LDP with rate function J, where 

I{a) = t{g t \a). 

Hence, under this superlinear time-scaling of the background Markov chain, the 
LDP for L M n [t ) is governed by a Poisson rate function with parameter g t . 

Example 5.2. Consider a modulated infinite-server queue (J, A, k , p), where the 
background process J is a standard Brownian motion W on [0, 00 ). By W we 
denote its restriction to the interval [0, i]. The sample paths of W are elements 
of Co[ 0 ,i], the space of continuous functions /: [ 0 , i] —> R with /( 0 ) = 0 . 

Equip Co[0, t] with the supremum metric. Of course, we may view the func¬ 
tion (f>t as a map from Co[0,t] to [0,oo) and this map is continuous under the 
supremum metric. 
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Scale A i —> nX and J K > J n , where J n is given by a linear time-scaling: 
J n {s) = W(s/n) for s > 0. Under this scaling, the arrivals are sped up linearly, 
whereas the time scale of the Brownian motion is slowed down linearly. 

This scaling resembles the scaling featured in 0. There, the authors consid¬ 
ered a modulated infinite-server queue under a linear scaling of both the arrival 
rate and the time scale of an irreducible Markov chain. The rate function ob¬ 
tained in [9] is given as the solution of a variational problem. We will obtain a 
similar result in this example. 

Since W is a Brownian motion, we have 

Schilder’s Theorem (cf. pTO] Th. 5.2.3]) states that ^=W satisfies an LDP in 
Co[ 0 ,f] with good rate function 

£(/) = A/o l/( s )| 2ds / S JTi([0,t]); 

loo else. 

Here, -ffi([0, t]) denotes the set of all absolutely continuous functions / € Co[0, t] 
that have square integrable derivative /. 

The contraction principle (cf. ]10j Th. 4.2.1]) now implies that (f>t(Jn ) satisfies 
an LDP with good rate function %Jj, where is given by 

= inf {£(/)l/ G #i([M)> Mf) = «}■ 

It follows from Theorem 13.41 that —M n (t) satisfies an LDP with rate function 
J, where / is given by 


1(a) 


inf 

767 ^ 0 ) 


[^( 7 ; a) + ^( 7 )]- 


Now recall that {?/> < 00} C 7 Z(t). Also observe that < 00} = Hi([0,t]) and 
that {i/j < 00} = { 0 t(/)|/ £ i?i([0,t])}. Then we may rewrite I as 

1(a) = inf [%;a) + ' 0 ( 7 )] 

l€{ip<oo} 

= . inf .MMf)>a) + ip(<j>t(f))]- 

Hence, / is given as the solution of a variational problem. 


6 Discussion and concluding remarks 

In this paper, we studied an infinite-server queue in a random environment 
and proved a full LDP for the transient number of jobs in the system. The 
proof of this LDP has two essential ingredients, namely the result that the 
transient number of jobs in the system has a Poisson distribution with a random 


18 


parameter and the assumption that the random parameter satisfies an LDP. 
Hence, the large deviations behavior of the random parameter seems to be the 
crucial factor that determines the large deviations behavior of the number of 
jobs in the system. 

The rate function corresponding to the LDP for the number of jobs is rather 
abstract. Nevertheless, we showed in the examples how to compute the rate 
function in certain specific cases. In particular, we recovered earlier obtained 
results for Markov-modulated infinite-server queues and strengthened these to 
a full LDP. Additionally, we proved LDPs when the background process has an 
uncountable state space. In all examples, knowledge about the behavior of the 
background process could be exploited to describe the rate function. 

There are several interesting topics for future research on the modulated 
infinite-server queue presented here. In this paper, we only looked at large 
deviations of the number of jobs at a fixed time t > 0. However, for certain 
applications it may be desirable to know the deviations over the whole time 
interval [0,i]. Therefore, it would be interesting to consider sample path large 
deviations. Also moderate deviations could be worth investigating, so as to 
bridge the gap between the central limit theorems and the large deviations 
results for modulated infinite-server queues. 

Furthermore, it could be very interesting to study other models. In par¬ 
ticular, more general arrival processes or service time distributions could be 
considered. As an example, it seems that the setup of [5] could be generalized 
to include a background process with countable state space and general service 
time distributions. Regarding general service time distributions, it should be 
mentioned that there might be some measurability issues when the background 
process has an uncountable state space. Finally, it would be interesting to see 
whether the large deviations results for modulated infinite-server queues carry 
over to modulated Ornstein-Uhlenbeck processes. To the best of our knowledge, 
this has not been investigated so far. 

Acknowledgement. This research has been partly funded by the Inter univer¬ 
sity Attraction Poles Programme initiated by the Belgian Science Policy Office. 

A Transient number of jobs in the system 

In this section, we provide the precise mathematical description of the model and 
determine the distribution of the number of jobs in the system at time t > 0, 
which is denoted by M ( t ). We mentioned in Section [1] that the steady-state 
distribution of the number of jobs in the system has already been determined 
for specific background processes in Model I and Model II. However, in this case 
we would like to determine the transient distribution given a general background 
process for the model described below, which generalizes Model I and Model II. 
Fortunately, the setup of our model is quite convenient and we may obtain the 
transient distribution without too much effort. 

By D([0, oo);£) we denote the space of cadlag functions from [0, oo) to £, 
where £ is a metric space with metric p. We define, in the usual way, a metric 
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d° on Z7([0, oo);£) that generates the Skorokhod J\ topology. (For more details, 
see Section II and references there.) 

Let (17, J- , P) be a probability space on which we have defined an indepen¬ 
dent, standard Poisson processes Y and an independent, cadlag stochastic pro¬ 
cess J with state space £. Assume that we have defined a collection of inde¬ 
pendent standard exponential random variables 2Ti, Z 2 ,... on this probability 
space. 

To modulate the infinite-server queue, we take continuous functions A: £ —> 
[0,oo), £ —> [0,oo) and k: £ —> [0,oo). More precisely, A (J) modulates 

the arrival rate, k(J) modulates the service requirement distribution and /i(J) 
modulates the server work rate. 

We define the modulated Poisson process Y via 



The process Y will be the arrival process. We denote the jump times of Y 
by 7i,72,... and the jump times of Y by ri,r 2 ,.... For convenience, we set 
7 q = tq = 0. The jump times T k and T k are related via T k = A ~(T k ) and 
Tk = A(rfc), where A (t) = fg A ( J(s)) ds and A _ (r) = inf{7 > 0 | A (t) > r}. 

Define the interarrival times ak = Tk — Tk -1 and 7Tk = Tk — fk —1 for k £ N. 
For later use, we note that 171 ,( 72 ,... is a sequence of i.i.d. random variables 
with a standard exponential distribution. 

At time t = 0 there are no jobs in the system. At each jump time of Y exactly 
one job arrives. Hence, the number of jobs that have entered the system during 
the time interval [0, t] is given by the (a.s. finite) random variable 1 {T- fe <*}- 

When job k enters the system at time Tk , it draws an independent service 
requirement from an exponential distribution with parameter n(J(Tk)) > 0, i.e., 
the service requirement of job k is given by Zk , where 



Job k leaves the system when its service requirement has been processed by the 
server, whose work rate is modulated by the background process J and is equal 
to n(J(s)) for s > 0. 

Hence, job k has both entered and left the system before time t > 0 if and 
only if Tk < t and Zk < /[r fc , fi(J(r)) dr. We get 



Note that M(t ) is a cadlag stochastic process. Because each Z k is strictly 
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positive with probability 1, it follows that 


M ^ 1 {z k <f t tATk ^J(r)) dr}) 

X-/( ^{ T k<t} ^{Zfc<K(J('Tfc))/ t t A M(J(>-))dr} )' 
fe= 1 v k k j / 


If J is deterministic, then it is relatively easy to determine the distribution of 
M(t). For instance, one may compute the characteristic function of M(t) via 
the following steps. 

Suppose that J(uj,t) = f(t) for all w € H and t > 0 for some function 
f € D([0, oo); £). For fixed n, g, / and t we define the functions g and h via 

g (s) = k(/(s)) f n(f(r)) dr, h(s) = 1 + [exp(i6>) - 1] exp(- ff (s)). 

J tAs 

Now we may write the characteristic function of M[t) as 


Eexp(i#M(t)) = Eexp^i0^^^1 {r(t < t} ~ ^-{z k<K (f( Tk )) (u ( /(r))dr}^ = 

oo / / n 

= E1 {-1>*1 + 5] E1 {r„<t;r„ +1 > t } exp id 71 - ^ t {z k <g(r k )} 
n— 1 \ \ k —1 



Clearly, El{ Tl>t } = e /o A (/(s))ds _ e A(t)_ yyg are computing the 

infinite sum above. Fix n € N and note that 


E1 {r„<t;r n+1 >t} exp ( id ( 71 - ^ 1 j^. 


k= 1 


{ Z k <g(T k )} 


— E l{r„<t;r„ + i>t} E 


exp 


(iO\n-Yt^ 


fc =1 


{ Z k <g(r k )} 


ri,r 2 ,... 


= E1 


71 

{r„<t;r„ +1 >t} ( 1 + [ eX P( i(9 ) “ l]e _!?(Tfe) ), 


fc=1 


because y and Zi, Z 2 ,... are independent. 

Next, observe that 

n n 

E1 {r„<i;r„ +1 >t} LI ( 1 + ~ l]e _9( ^ fc) ) = El {rn < t;Tn+1>i} h(r fc ) 

k =1 fc=l 
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and 


E1 {r„<t;r n+1 >t} fj h{r k ) = if t{ Tn <t} ( f[ K T k) ) E i 1 {cr n+1 >t-r n } \ Tl,...,T n ] 

fc=i V \fe=i / 


= E ( l{ Tn <t} ( h{T k ) 


0 -(A(t)-A(r„)) 


\k-l 


For convenience we write x k = x\ + ■ ■ ■ + x k . We have 

k ( %»<*> f n m^)V a(tb)> ) = 


Kk=l 


— E l{T„<A(t)} 


n^ A ( r fc)) 


\k— 1 


r A(t) r A(t)-x+ r A(t)-x+_ 1 


' £1—0 j £2—0 


• Xn— 0 


IIMA- (x k )) da;„ ... dxi 


fc=1 


fA (t) r A(t) f A(t) n 

= / / ••• / n^ A (Vk))dy n ...dyi 

Jy i =0 Jy2=yi Jy n =y n -1 k=i 

pt pt pt n 

= / ••• / JJ [h(zk)X(f(zk))] dz n ... dz\. 

J zi=0 J Z2—Z\ J Z n =Zn-l 

Now note that for an integrable function g we have 


g(s) ds 


uo 


= n! 


! Z\—0 j Z 2 =Z\ J Z n —Zn- 1 


]~[ g(z k )dz n ...dzi. 


k =1 


As a result, it holds that 

rt rt rt 


U [h(z k )X(f(z k ))] dz n ... dzx = 

J Z1=0 j Z2 — Z1 j Zn = Zn -1 

= ~\ [ h(s)X(f(s))ds 
U o 

n 11 / /** 

= Y, fcI A ^) fc (n _ fc)! (j ex PW - M J q Kf(s))e~ 9{s) ds 


i — k 
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Now we may write 


Eexp(i 6M(t)) = 

oo / / n 

= ®1{ti> 4} + ®^{Tr.<t;r n+ i>t} exp ( ( n — ^ ^{Z k <g(T k )} 


n—1 


V 


k=l 


= e~ A ^ 


OO n 1 1 

E«- A(,, Es A w‘ 


n—1 
oo n 


k—0 


\e ld - l] J X(f(s))e s(s) ds 


(n — k)\ 

OO Tl 1 , / /.t \ ’ 

= e_A< ” ? ? EF A(t >‘(^iy 0” - ') l A ‘) 


n—k 


i-k 


n—0k—0 
oo oo 


= e- A « 


^ ^ 11 / r 1 \ 

EEy A W^ ([e W -!] / A(/( S ))e-^)d S ) 

fe = 0 ra =0 A J ° A 

= exp ^[exp(i0) — 1] J A(/(s))e _K ^ s ^ ^(-f( r )) dr ds ^. 


Hence, in this case M(t) has a Poisson distribution with parameter where 

Mf) = f 0 f Mf(s))e- K(f( * ))f > (f(r))dr ds. 

Now suppose that J is not deterministic. Then J is a random element of 
D([0,oo );£). In this case, we may use the independence of J and standard 
arguments to obtain that 

Eexp(i 6M(t)) = EE[exp(i0M(t))| _F d ] 

= Eexp^[exp(i0) - 1 ]J X(J(s))e~ K(J(s)) J> (J(r))dr ds^j. 


We summarize our findings in the following lemma. 

Lemma A.l. Under the stated conditions, M(t) has a Poisson distribution 
with random parameter 


pt 

MJ)= / A(J(s))e- K(J(s)) ^> (J(r))dr ds. 

Jo 

Consequently, if we scale X(x) i->- nX(x) and J K > J„, then the number of 
jobs in the system M n (t ) has a Poisson distribution with random parameter 
n<f>t{J n )- This observation is crucial for the proof of the LDP for -f M n (t). 

B Continuity in Skorokhod space 

Let £ be a metric space with metric p. Let D([0,oo);£) denote the space of 
cadlag functions /: [0,oo) —» £, i.e., lim s j. t f(s) = f(t ) and lim s ^ t f(s) exists in 
£ for every t > 0, where lim s ^o f(s) := /(0) by convention. 
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Define a metric d° on D([0,oo);f) via 


d°(f,g) = inf 

AgA 


7 (A) V / e U d(f 1 g 1 \,u)du 


Here, A denotes the space of increasing homeomorphisms of [0, oo), 


7 (A) = sup |log(A(f) - A(s)) - log(t - s)| 

£>s>0 


and 


d(f,g,\,u) = sup [1 A p(f(tAu),g(\(t)Au))}. 

££ [0,oo) 

The metric d° induces the Skorokhod Ji topology. For more details, see m or 

m- 

Definition B.l. A function f c £ D([0, oo);£) is called a piecewise constant 
function or a step function if there exist n £ N, finitely many time points 
0 = to < 1 1 < ... < t n < oo and ao ,...,a n £ £ such that f c (t) = ai for 
t £ [tj, ti+ 1 ) and i = 0 ,..., n — 1 and f c (t) = a n for t £ [t n , oo). 

The set of step functions in D([0,oo);£) is denoted by S([0,oo);£). 

Proposition B.2. Let f £ D([0, oo); £). For all T > 0 and e > 0 there exists 
a step function f c £ <S([0, oo); £) such that 

sup p(f(t),f c (t)) <e. 
te[o,T] 

Proof. This is derived in the same way as [16[ Th. 12.2.2]. □ 

Corollary B.3. The set <S([0, oo); £) is dense in D{[ 0,oo);£). 

Consequently, every continuous function on D([0,oo);f) is completely de¬ 
termined by its behavior on the set of step functions. 

Let A: £ —> [0,oo), k: £ [0,oo) and p: £ —> [0,oo) be continuous. For 

t > 0, we would like to show that the function <j) t : D([0, oo); £) —> [0, oo) defined 
via 


pt 

Mf)= / A(/(s))e- K(/(s)) ^> (/(r))dr ds (14) 

Jo 

is a continuous function. 

First, we observe that the map c*: D([0,oo);£) —> D([0,oo);R) defined 
via C\(f)(t) = A(/(f)) is continuous, because A is continuous. Similarly, the 
functions c K and c p are continuous. 

Next, let f,g £ D([0,oo);R). Then pointwise multiplication of / and g , 
defined via ( fg)(t ) = f(t)g(t). This is a measurable map which is continuous 
(/,<?) if / or g is continuous (cf. I FT Th. 4.2]). 
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Finally, let / £ P([0, oo); R). Then the map ip: D([0, oo); R) —» P([0, oo); R) 
defined via ip(t) = J 0 * f(s) ds is continuous. This follows almost immediately 
from the definition of ip and the characterization in in Pr. 3.5.3]. 

Now note that the sequence of functions {A(/ n )} n6N is bounded in the sup 
norm over [0, £] if /„—>•/ in D([0, oo); £). Hence, it suffices to show that 

f e -K(/„(s))/ s V(/n(r))dr f e ~K(f(s)) f* fi(f(r))dr fa 

JO Jo 

as /„ —> f in P([0, oo); £). But this follows from repeated applications of the 
first three observations. 

Hence, the map (pt must be continuous. Note that continuity of A, n and y 
is crucial to obtain this result. We summarize these findings in the following 
lemma. 

Lemma B.4. Let A: £ —> [0,oo), k: £ [0,oo) and /x: £ —> [0,oo) be contin¬ 
uous. Then the function <p t : P([0,oo);£) —> [0, oo) as defined in equation (1141) 
is continuous. 


C Properties of Poisson random variables 


For 7 > 0, let Po(y), Pi ( 7 ), P 2 W, ■ • ■ denote a sequence of i.i.d. random variables 
that have a Poisson distribution with parameter 7. In this section, we will fix 
an arbitrary i£l, <5>0, A>0 and e > 0 and define A)T = max{0, A — e} and 
A+ = A + e. Recall that B+(A, e) = B( A, e) D R + . 

We would like to prove a large deviations lower bound for 

1 (1 " 

liminf inf — logP — Pi( 7 ) S B(x, S) 

n -¥00 7 es + ( a,c) n \ n ' 

\ i—1 

Of course, the difficulty here is the presence of the infimum over a range of 
parameters. We will show in Proposition lC.il that this infimum may be taken 
over certain restricted subsets of P+(A, e). For each of these subsets we will 
provide a large deviations lower bound, from which we will derive a lower bound 
when the infimum is taken over B + ( A, e). This is the content of Proposition lC.31 

Proposition C.l. For all x £ R, <5 > 0, A > 0 and e > 0 it holds that 


inf 

7 eB+(A,e) 


1 n 

~^2 P i{l) e B(x,S) 
n 


inf 

7 e(B(A,e)nB[a:,5])u{Ar,A+} 


1 " 

-E P *(7 )£B{x,6) 


Proof. Let 0 < 7_ < 7+ < 00. For y £ R it holds that 


P(P 0 (7 + ) =y)> P(P 0 (7-) =y) if y> 7+ > 7 - ( 15 ) 
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and 


P(P 0 (7+) = y) < P(P 0 (7-) = 2/) if 7+ >7- >V- (16) 


Because we are working with i.i.d. Poisson random variables, we may write 

P^y^-Pj( 7 ) € = P(P 0 (« 7 ) e (n(x - 5),n(x + 6))). (17) 

Now the statement of the proposition is an easy consequence of the equations 
(HU), USD and (ED combined. □ 

Proposition C.2. Let x £ R and 5 > 0. If B + (x,S) ^ 0, then 

lim inf — logP ( — Pi( 7 ) £ B(x, 5) ] =0. 
n^oo-yeB + [x,5] U \ n ~~[ ) 

Proof. For a Borel set A C R, define p n (A\^) = P(1 Piil) £ 7l). Now 
suppose that B + (x,8) / 0. Then the diameter of B + (x,S) is strictly positive 
and bounded above by r = min{2i5, x + c)}. 

Let N r £ N be such that jj- < %. Then for all n > N r and 7 £ B + [x, 5] we 
define 7 “ = l|n 7 _l> In = £T n 7 l and 

7n = min {{7n >7« } r\B(x,S)}. 

Then max {|7 - 7 ~|, |7 - 7 +|} < 1 < § and p n (B(x, < 5 )| 7 ) > Pn({7n}l7) for 
each n £ N and each 7 £ B + [x, < 5 ]. Using that 77! < n n + l / 2 e ~ n + l j we get 


Pn((7n}|7) > 


77,7 


777+I 

> 11 - 




e (nj*) 


- 1/2 


1 


) n(x+S) 

e~ 2 (n(x + 5))~ 1/2 


n{x + 8) + 


for each 77 £ N and each 7 £ B+[x, <5]. This implies the statement. □ 

Combined with Cramer’s Theorem in R, the two previous propositions enable 
us to prove the following large deviations bound. Note that we prove an equality 
rather than an inequality and that the limit exists. 

Proposition C.3. For all x £ M, <5 > 0, A > 0 and e > 0 it holds that 


lim inf — log P 

n -¥00 7 gB + (A,e) 77 



Z p ‘(7 )eB(x,8) 


min 

7<E{A e -,A +} 


a 


inf 

EB(x,S) 


t(r,a) 


(18) 
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Proof. Define p n (A | 7 ) = P(i Y^i =1 -^( 7 ) £ -4) for Borel sets 4cl and C = 
(B( A, e) fl B[x, <5]) U {A 7 , A+ }. Thanks to Proposition [Cj] we may write 

lim inf — logp n (B(x, S) | 7 ) = lim inf — logp„(-B(o;, 6) I 7 ). 

n->oo 7 gB + (A,e) 77 n->oo 7 gC 77 

It follows from Proposition IC.2I that we may restrict the infimum to the set 
{Ar, A+}, so 

lim inf — logp„(_B(:r, S) | 7 ) = lim min — logp n (B(x, 5) | 7 ) 

n—too 7 gC 77 Ti-too 7 g{A7 ,A+} ?7 

= min lim — logp n (B(x, 5) \ 7 ) 
7e{Ar,A+} n ^°° n 


= min — inf £( 7 ; a) . 

7 e{Ar,A+} L a&B(x,6) 

The last equality is an application of Cramer’s Theorem for i.i.d. Poisson random 
variables; the limit exists because B(x,5) is a continuity set for the Fenchel- 
Legendre transform corresponding to a Poisson distribution. □ 

As shown in the inequalities (0 and m, the Fenchel-Legendre transforms 
corresponding to Poisson distributions are nicely ordered in some sense. This 
property leads to the following propositions. Their proofs are elementary but 
tedious and are therefore omitted. 

Proposition C.4. Let F Cl be closed and define f: [0,oo) —> [— 00 ,0] via 

f(l) = - inf %;a). 

If F C (—oo,0), then f = — 00 . If F fl [0, 00 ) ^ 0, then f is real-valued and 
continuous on (0, 00 ). Additionally, lim 7 j,o /(y) = /(0), where /(0) = 0 if 0 £ F 
and /(0) = ooif0^F. In any case , / -1 ([a, b\) is closed for all a,b £ (— 00 , 0] 
with a < b. 

Proposition C.5. Let IZ C [0, 00 ) be a non-empty, closed set. Let ip: R —> 
[0, 00 ] be a lower semi-continuous function. Then the function I: K —> [0,oo] 
defined via 


1(a) 


inf [£( 7 ; a) 

7G7c 


^(7)] 


is a lower semi-continuous function. 
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